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Abstract. We present numerical simulation studies of 2D reduced MHD equations 
investigating the impact of the electronic /? parameter and of curvature effects on the 
nonlinear evolution of drift tearing islands. We observe a bifurcation phenomenon 
that leads to an amplification of the pressure energy, the generation of E x B poloidal 
flow and a nonlinear diamagnetic drift that affects the rotation of the magnetic island. 
These dynamical modifications arise due to quasi linear effects that generate a zonal 
flow at the onset point of the bifurcation. Our simulations show that the transition 
point is influenced by the /3 parameter such that the pressure gradient through a cur- 
vature effect strongly stabilizes the transition. Regarding the modified rotation of the 
island, a model for the frequency is derived in order to study its origin and the effect 
of the /3 parameter. It appears that after the transition, an E x B poloidal flow as well 
as a nonlinear diamagnetic drift are generated due to an ampliflcation of the stresses 
by pressure effects. 
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1. Introduction 

In tokamak and space plasmas, confinement can be affected by instabilities and 
in particular, at resonant surfaces, magnetohydrodynamics activity can lead to the 
generation of magnetic islands reaching a macroscopic width. Solar flares jT], energy 
release events in the geotail [2] or tokamak internal disruptions, also known as sawtooth 
oscillations, are linked to such reconnection phenomena. Diamagnetic effects and self- 
generated zonal flows can modify the saturated island width via bifurcation mechanisms 
[2]. The rotation frequency of the island can also be nonlinearly affected with a strong 
dependence on the transport coefficients and on the competition between the Reynolds 
and Maxwell stresses [1]. This can have a significant physical consequence, for example 
in a tokamak, where such a nonlinear effect on the rotation can lead to a slowing down 
of the plasma through locking to the resistive wall producing in turn a degradation 
of the plasma and/or triggering a transport barrier [5]. Likewise curvature effects can 
also modify the nature of island dynamics. Magnetic islands can in particular coexist 
with pressure driven intabilities such as interchange modes and/or turbulence. Several 
experiments report the coexistence of turbulence and MHD activities showing some 
correlated effects [6l |7]. Numerical studies of the interaction between double tearing 
modes and micro-turbulence to delineate the interaction between zonal flows and the 
latter in the growing phase of the double tearing instability have also been performed 
in [H]. More recently, in [H] an investigation of the interaction of a 2D electrostatic 
turbulence with an island whose dynamics is governed by a generalized Rutherford 
equation has been carried out. However the study neglects the potentially stabilizing 
influence of the magnetic structure on the turbulence precluding thereby any multi- 
scale interaction between MHD and turbulence. In this paper we study the dynamics 
of a magnetic islands in the presence of interchange effects but limit ourselves to the 
situation where the interchange modes in the system are linearly stable. We find that 
the whole system does not generate turbulence in the nonlinear stage but exhibits a 
complex dynamics arising mainly due to quasilinear effects. Our investigations are 
based on linear and nonlinear simulations of a set of reduced fluid equations (a three 
field model) through which we examine the origin and the influence of zonal flows on 
the magnetic island dynamics in the presence of interchange effects. The magnitude of 
the pressure gradient appears to be a key parameter of the dynamics controlling both, 
the generation of the zonal flow and the development of a nonlinear transition in the 
system. 

The paper is organized as follows. In section II, the model equations are introduced. 
In section HI, a linear analysis of the model is done in order to understand the role 
of the equilibrium magnetic field in the stabilization of the electromagnetic interchange 
modes. In section IV, the description and the analysis of the dynamics are done. In 
section V, the origin of the island poloidal rotation is investigated. Section VI presents 
a summary and conclusions of the paper. 
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2. Model system 

Our model system is a three fields model corresponding to a reduced magnetohydrody- 
namic description of the fiuid equations [TU] and which provides a minimal framework 
for including both the interchange and the tearing mode phenomena in a plasma. The 
model consists of a set of three coupled equations for the electrostatic potential 0, the 
pressure of the electron p and the magnetic fiux ip. We suppose that the magnetic field 
is dominated by a constant component Bqz along the z-direction. The time evolution of 
the three fields are described by: 

dtyl<j) + [0, Vi</.] = VltP] - K.dyP + z/vl</), (1) 

dtp + = - - K2)dy(l) + K2dyp^ + [tp, ViV] + X±Vip, (2) 

dtijj + [(f) - p,ip] = -v^dytp + rjV^ip, (3) 

where v^, = 2nt^L ■ "^^^ sum. of the electron and ion momentum evolution equations 
leads to the plasma equation of motion, Eq. ([T]), where v is the viscosity. Eq. ([2]) comes 
from the energy conservation equation where xi. is the diffusivity. Eq. ([3]) is Ohm's 
law (electron parallel momentum equation) with rj being the resistivity. /3 = ^^^"^^^ is 
the ratio of the electron thermal energy to the magnetic energy (j>o being the amplitude 
of the equilibrium pressure), Lp is the pressure gradient length, is a magnetic shear 
length, Rq is the major plasma radius, Vti = is the ion cyclotron frequency, and Ta 
is the Alfven time. Equations ([UlS]) are normalized as follow: 

>t, (4) 

^^A, 7 5 ><P, -j-^P^P^ (5) 



L±Boz L±vaBoz L^po 

where va = Bqz/ fionrrii = L±/Ta is the characteristic Alfven speed. parameters are 
linked to the curvature and to the pressure gradient (^i = 2QiTA^ and K2 ^ ^^^"'^ 



so these parameters control the interchange instability. On the other hand, in Eq. ([2]), 
the tearing mode dynamics is controlled by the coupling parameter = gj^rr- More 
precisely, this parameter controls the coupling between pressure and the magnetic fiux. 
The nature of the linear and nonlinear dynamics of the magnetic island depends strongly 
on the strength of the coupling. For a high /3 plasma, since the coupling is strong, the 
pressure and the magnetic fiux control the island dynamics, whereas for a low /3 plasma, 
the island dynamics is governed by the interaction between the fiow and the magnetic 
fiux. In our model we assume the electron temperature to be constant and the ions to 
be cold. The cold ion limit is physically realistic since the ion temperature does not 
significantly affect the stability of the tearing mode. As a further simplification we have 
also neglected the parallel ion dynamics in the energy balance equation Eq. ([2]). Eqs. ([1]- 
[3]) are solved numerically using a finite difference scheme in the x direction, including an 
Arakawa algorithm [11] for an accurate conservation of the Poisson brackets [., .] and a 
pseudo-spectral method in the y direction, including an appropriate de-aliasing scheme. 
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C2 


0.001 


2* 10-3 


3.3 * 10-3 


0.005 


10-2 


1.67* 10-2 


0.015 


3* 10-2 


5 * 10-2 


0.025 


5* 10-2 


8.33 * 10-2 



Table 1. Effect of the (3 parameter on v* and C. 



3. Nature of the tearing modes and the influence of the curvature 
parameter 

We now study the influence of the interchange mechanism on the magnetic reconnection 
when the gradient scale length of the pressure, Lp, is of the order of the size of the island: 
we set Lp = L±. We are interested in large islands, i.e islands with widths w such that 
a > w ^ Ps where a is the minor radius and ps is the hybrid Larmor radius (p^ = Cg/fli, 
where Cg is the ion sound velocity), and we have chosen L±_ = 0.24m. The numerical 
values of other parameters are taken to be Rq = 2.24m and ^2^X4=0.5. These numerical 
values are typical of the TORE SUPRA device for an island width of about 1/3 the 
minor radius and lead to Hi ~ 0.11 and K2 ~ 0.36. The widths of the numerical 
integration box are set to = 2ttL± and Ly = 5'n'L±. The values of the coeflicients 
C{f3) and = kyV^, are determined for four different values of /3 in the range of 10-^ 

to 2.5 X 10"2 (see Table ([1])). Transport coeflicients (i/, 77,%) are all set to 10-^ which 
correspond to renormalized coeflicients to include eflects of microscopic turbulence [12] . 
The equilibrium magnetic fleld Boy = Boy, based on the Harris current sheet model |13| . 
is chosen to be of the form, 

Bo{x) = tanh f^LlllIl\ . (6) 



at 

The parameter at = 0.75 controls the width of the proflle, iPq{x) = Bo{x). With such a 
profile, the parameter A' (the tearing mode stability index) can be explicitly computed 
taking into account the boundary conditions, for modes evolving slowly on the Alfven 
time scale and further neglecting the viscous and interchange corrections. Introducing 
k = atky, we have 



atA' = 2(1/ k-k) (7) 



where 



/Lx/i2at) 
rfyexp(2%)/(l + tanh(y)/A;)2 (8) 



is a correction linked to the finite radial distance of the walls. 



Let us investigate the stability of the modes modeled by eqs. ([IHS]) with the given 
numerical values of the parameters in the presence of such an equilibrium. Figure ([T]), 
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Figure 1. Linear growth rate 7 [r^ ] versus poloidal mode number ky for a simulation 
with /3 = 10-=^. [Left] w* 2.10"^, C = 3.33. 10"^, ni = 0.1071, ^2 = 0.3571 
and ^ X± = V = 10-"^. [Right] = 10"^ C = 2.10^3, = 5, K2 = 0.3571, 
fi = X± = 10"^ and = 10""^. 



shows the growth rate of the electromagnetic interchange and the tearing branches, as 
functions of the poloidal wave number, for the parameter values given above. The left 
graph has been obtained using the parameters chosen in this work with /3 = 0.001. 
Tearing instability has the largest growth rate at ky""" = 2t[/LY = 0.4, for which 
A' = 7 and '-ftcn.v ~ 0.0042. This is clearly smaller than the one we would obtain 
in the classical tearing limit , i.e if we would have set all the parameters to zero 
except rj (7tearciass = 0.0072). It is instructive to note also that in this parametric 
regime the interchange branch is stable for any wave number. From an electrostatic 
point of view, with such parameters, interchange would have been unstable for ky < 8 
and would have given a scale separation between both instabilities {ky'"'^'"'/ky'"" ~ 7, 
^int oiec^yoar ^ -^y j^g^ fg^^^g ^j^g tcarlug brauch. Linearisation of eqs. ([T]-[2]) in the 
vicinity of the resonnance shows that curvature effects weakly modify the growth rate 
if KiVi,{ky/kx)'^ /'^'^ <^ 1 and k,2 ^ 1, which is true in our case. Considering the fact 
that the linear regime is also not controlled by viscous phenomena jT^, it follows that, 
linearly, this system develops approximately drift tearing modes. The actual nature of 
these modes is controlled mainly by the ratio P/Pqi and 7/0;^, where P = v/rj = 1 is 
the Prandtl number and Per = {^'{ri/kyY^^y^^. In our cases, for any (3, the first ratio is 
always smaller than 1. It implies that when > 1, one gets the visco-tearing regime 
with a growth rate scaling law 7^,^ ~ 0A7A'ri'^^^P~^^^kl^^. When 7/0;^ < 1, we recover 
the visco-drift-tearing regime with the growth rate 'y^dt |15| . For instance, for /3 = 0.001 
and ky = 0.4, we have 7^,4/0;^ ~ 6.5 > 1 and 7 = 0.00042 ~ 7^^ = 0.00052. For ky = 1.2, 
we have -f^t/uj^ ~ 0.2 < 1 and 7 = 0.00035 ~ -f^dt = 0.0004. 

The right graph of Figure ([T]) shows that there exist regimes where the interchange 
branch is unstable and has the largest growth rate at small scales. The study of such 
regimes is out of the scope of this paper. We also remark that the instability does not 
necessarily develop in the vicinity of a resonant surface, but in that case, the effect of the 
magnetic field on the stability of interchange like modes can be investigated by setting 
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Figure 2. Stability of electromagnetic interchange modes versus -0q, far from the 
resonant surface. Same parameter values as figure 1. 
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Figure 3. Time Evolution of the Energies for /3 = 10 



ipQ = Bq = Cte and using some Fourier analysis. Figure (|2]) shows the hnear growth rate 
of interchange modes versus ip'^. As is well-known |16| . the equilibrium magnetic field 
stabilises the interchange modes and, in our case i?o = 1, this is clearly stable. We can 
therefore expect that in the initial phase the growth of the magnetic island is weakly 
infiuenced by interchange parameters. 



4. Nonlinear generation of a strong zonal flow 

4.1. Description of the nonlinear evolution of the system 

To characterize how the pressure gradient affects the evolution of a magnetic island, 
linear and nonlinear self consistent numerical simulations have been performed. A grid 
number of = 128 is chosen for the radial direction and Uy = 128 for the poloidal 



Figure 4. Snapshots of the fields -0, and p before transition for /3 — 0.001, at 
t = SOOOta. 



direction (equivalent to 48 modes in this direction, including dealiasing). The energy 
conservation relation derived from Eqs. ([IHS]) is 

^ {E^ + Ep + Ek) = -ri< f >-u< > < |Vp|' > +S , (9) 

where E^ = 0.5 <\V{ij - ipo)f >, Ep = 0.5 <p^> /C^ and Ek = 0.5 <|V0|^> are 
respectively the magnetic energy, the pressure energy and the kinetic energy of the 
fluctuations. The brackets < . > mean here an average over the simulation domain. S 
is the source term linked to the curvature and the pressure gradient, proportional to 
the radial pressure flux , S = —as < pdycf) > with = ^(1 — K2) + > because 
< 1. Note that a local flattening of the pressure by radial exchange of pressure cells, 
gives a fluctuation 6S < 0. Moreover, the interchange source term S is not modified 
by the generation of zonal flow. Figure shows the time evolution of E^, Ep and 
Ek for the parameters chosen in this work with /3 = 10~^ as well as the corresponding 
E^ and E^'^, for a classical tearing mode (i.e p = and Ki = 0). In comparison with 
the evolution of a classical tearing mode, four regimes are observed in the nonlinear 
simulations of a magnetic island in the presence of the interchange term. First, there 
is a linear regime where the magnetic island is formed. Second, the system reaches 
a quasi-plateau phase. Then, a transition occurs and as it will be shown later this 
is linked to the interchange parameters. Finally, the system reaches a new saturated 
state. During the first two phases, the evolution of the energies is not strongly affected 
by the presence of the curvature terms. The evolution of the magnetic island follows 
closely the time trajectory of an island driven by a tearing instability. However, at 
t^, = 13200 ta, a transition occurs. Figures (jl]) and ([5]) show snapshots of the fields ip, 
and p respectively before and after the transition. The two dimensional profiles of the 
pressure and the electrostatic potential (represented through isocontours) are strongly 
affected by this transition. After this phase, the structure of the mode changes and a 
flattening of the pressure is obtained. 
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Figure 5. Snapshots of the fields and p after transition for /3 = 0.001, at 
t = 18000rA. 
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Figure 6. Spectra before (a) and during the transition (b), at respectively t = 12000rA 
and t = 14500t_a. 

4 .2. Origin of the transition 

In order to understand the origin of the transition and to characterize the structure of 
the electrostatic potential after the transition, we first assess the importance of small 
scales. The spectra before the transition (a) and during the transition (b) are shown in 
Figure ([6]). An equipartition between the energy of the magnetic flux and the energy of 
the pressure is observed at large scales 0.8 < ky < 5 whereas there is an equipartition 
between the energy of the pressure and the kinetic energy at small scales 8 > ky > 14. 
We observe that these properties continue to persist when the transition occurs, the 
energy level of the modes ky > 1.2 being roughly unaltered. We also note that as 
apparent in Figure ([6]) only the large scales are affected at the transition. A detailed 
analysis shows that the transition occurs when the kinetic energy of the mode ky = 
becomes equal to the one of ky = 0.4 (mode 1). This suggests that the dynamics 
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Figure 7. Runs with same parameters as figure (3) : [Left] Time evolution of the 
energies with ki = K2 = (4 modes runs) . [Right] Time evolution of the kinetic energy 
of modes 0, 1 and 2 along with the corresponding classical tearing case. 



of the structure of the electrostatic potential is quasi-linear and that the transition 
occurs when the mode ky = becomes energetically dominant. To delineate the quasi- 
linear nature of the magnetic island dynamics, we have performed a simulation with 
only four poloidal modes and with the same parameters (/3 = 10~^). Figure ([7]) shows 
the time evolution of the energies for this simulation. The comparison with Figure 
(12]) demontrates that one needs only four modes to describe the time evolution of the 
energies. Therefore the dynamics of the system is the result of quasi linear effects. 
Further, the time evolutions of the kinetic energies of modes 0, 1 and 2 are presented in 
the right hand side panel of Figure ([7]) and compared with the evolution of modes in a 
classical tearing run, i.e without pressure effects {p = and Ki = 0). In the latter case, 
the growth is driven by the mode 1, the transition does not occur and the mode is not 
generated. Conversely, in the case where the pressure effects are included, the transition 
occurs and the mode ky = 0, i.e zonal flow, is strongly generated. It is also the first to 
be amplified exponentially, at the beginning of the transition. This suggests that the 
transition is linked to a strong amplification of the zonal flow. Nevertheless, when we 
perform the same run, suppressing artificially the mode 0, i. e the zonal flow, we find that 
the transition still occurs, roughly at the same time, but with a weaker amplitude. This 
suggests that althouh zonal flows play a predominant role they are not the sole factor 
responsible for the destabilizing mechanism. Indeed we can show from the analysis of 
the snapshots of 0, p and ip given by Figure ([H]) that the mode plays an important 
role in the triggering mechanism of the transition phase. We can observe that during 
the transition, the pressure cells are crossing the resonant surface at the current sheet in 
both directions, producing at the end a modification of the mode structure. Clearly the 
potential structure suggests that the crossing results from an advection by the flow. Let 
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Figure 8. Snapshots of the fields ■0, ip and p for /? = 0.001, at t = ISOOOr^i. 



US stress that during this phase, the classical tearing picture of an incoming flow from 
the sheet into the island is no more valid. Between t = ISOOOr^ and t = ITOOOr^, S{t) 
decreases because the reorganisation of the cell is radially equivalent to an exchange 
of pressure cells with the gradient of the pressure fluctuations being outward. From 
an energetic point of view, see eq. ([9]), Ki and K2 have a negligible effect because the 
dominant contribution in the interchange source term is linked to w^/C^. In Figure 
d?]), the results of a simulation with four modes and without the curvature terms are 
presented (ki = ^2 = 0). We flnd that the transition does not occur. At least, it does 
not occur at 2.5^^,, showing that a more complex mechanism due to curvature terms 
might be at play. We next investigate the origin of the zonal flow. 



4.3. Origin of the strong increase of the zonal flow 

We investigate the origin of the zonal flow that is generated whenever /3 7^ by 
considering separately the energy transfer from the tearing mode to the zonal flow 
and from also the other modes, in particular, the small scales. 

Following [S], the equation for the flow energy feeding the mode = m2'K / Ly can 
be written as, 

^E^. = Ti + T*^ + T^ + T^^ + T^' (10) 



dt 



where 



T« = -/dx^„(l*.-l)..(R«y.>olds stress contribution). 
r« = (Maxwell stress contribution). 

T^^ = J dx^trnilipo, j])m (line bending term contribution), 
dx(j)m{[ip, jo])m (kink term contribution). 
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Figure 9. (3 



10-3; Time evolution of , T^^ , T^^ and T^^ . 



Here 



, where 0^ is the m Fourier component of 



In order to understand the origin of the generation of the zonal flow, it is useful to 
project equation (fTOj) on the mode m = 0: 



M 



rpLB I rpKI 



The curvature term does not directly feed the zonal flow. In Figure (Q, the time 
evolutions of T^, Tq'^ , Tq^ and T^^ for the simulation with /3 = 10"'^ are presented. 
The contributions of the line bending and the kink terms are very weak. However, at 
the transition, there is a strong generation of the Reynolds stress and the Maxwell stress 
contributions. To proceed further, it is useful to separate the contributions from zonal 
flow (m = 0), the mode m = 1 and other modes for each of the transfer functions. 
Let us introduce 0>i = X]m>i ^m- We can define three contributions in each transfer 
function. For instance, in the case of Ti?, we have: 



where Tj^ 



and T^^ 



rpM 



>lm 

0. Let us 



= / dx(l)m{[(pO,(^o])m, T^rri = I ^X^m ( [01 , Wi] ) 

/ dx(f)m{[4>>i,uJ>i])m- Clearly, by definition T/J = and = 
focus on the energy transfer to zonal flow, neglecting the weak contributions of the line 
bending and the kink terms. Equation (fTTl) then becomes: 

j^Eo = TP, + T^,, + T^^ + T^{, (13) 

Using the above prescription we have checked that the main contibution to the Reynolds 
and Maxwell stresses comes from the mode m = 1 while the contributions of the small 
scales are weak. 



4-4- Effect of the /3 parameter on the nonlinear dynamics 

The left panel of Figure f llOp presents the transition time and the time where the 
first quasi plateau saturation occurs for various values of /3. We note that the time 
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Figure 10. Effect of the /3 parameter and of K2 on the transition : [Left] Time at the 
beginning of the first quasi plateau phase and time at the beginning of the transion 
versus /3. [Right] Time evolution of the kinetic energy for simulations with (3 = 10~^ 
but with different values of ^2- 



corresponding to the first quasi-plateau phase does not depend on the /3 parameter. 
This is in agreement with the results of Figure ([3]) which shows the time evolution of the 
energies for two cases namely the classical tearing case (/3 = 0) and the /3 = 10"'^ case. 
However the saturation time depends strongly on the /3 parameter. For a regime where 
the pressure effects are strong, i.e for a high value of /3, the transition occurs quickly 
whereas, for low /3 regimes, the transition occurs later. Further, in the right panel of 
Figure (fTOj) . the effect of the interchange parameter K2 on the transition is shown for a 
simulation with four modes and with /3 = 10~^. The transtion time depends clearly on 
the K2 parameter. It tends to stabilize the first plateau phase. 

To summarize, the nonlinear transition results from quasilinear effects. Zonal flow 
amplification at the transition is due to the energy transfer from the mode m = 1 to 
the mode m = through mainly the Maxwell stress. The curvature term K2 linked to 
the interchange effect does not directly feed the growth of the zonal flow. However, we 
have shown that this term controls the transition time t*. The transition leads to an 
effective radial exchange of pressure cells generating an outward mean pressure gradient 
of fluctuations. The shape of the pressure structure after the transition implies that a 
diamagnetic velocity a)^ has been nonlinearly generated, driving a rotation of the island. 
This driving is of course in competition with the zonal flow. Let us in the next section 
analyze quantitatively the island poloidal rotation. 
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5. Study of the island poloidal rotation 

5.1. Model for the island rotation frequency 

Following |1] where a study of the rotation frequency of the island has been done for the 
case of a drift-tearing mode, we investigate the origin of the magnetic island poloidal 
rotation. Let us project the Ohm's law (Eq|2]) on the mode m = 1 assuming that for 
the mode fc^, ipky{x,y,t) = ipky{x)e^'''^^e~^'^^, where the real part of 7 is the frequency 
of the island rotation and the imaginary part is the linear growth rate of the island. 
Neglecting the nonlinear contribution of the modes ky > ki, we obtain the expression 
for the rotation frequency u of the island 

U = U* + + UexB + L^o + Lr, (14) 



where 



uj* = kiv*, 

uj* = - kid^po, 

Coexb = kidx(j)Q, 



L^o = - Re (^kitpQ — 
= -Re I IT] 



{dl - kl) (x) 



^fci (x) 



are respectivly the linear diamagnetic drift, the nonlinear diamagnetic drift, the 
contribution of the equilibrium magnetic field and the contribution of the resistivity. In 
general, each term of Eq.( lT4|) is not a constant inside the current sheet 6, so we consider 
their radial averages over the current sheet to contribute to the rotation frequency. Eq. 
f fl^ becomes : 

< U >5 = < + U* >5 + < UExB >5 + < L^:o >5 + < Lr, >5 (15) 

where < . >s means an average over the current sheet. 

In Figure (fTTj) . the time evolution of the poloidal position is presented for a non- 
linear simulation with /3 = 0.025. The dynamics of the energies for a simulation with 
P = 0.025 has the same behaviour as the one obtained in Figure ([3]) with /3 = 10~^. 
However, as shown in the left panel of Figure f fTOj) . with such a high value of /3, the 
transition occurs earlier around t = 5600rA. In Figure (ITT]) a comparison with the is- 
land position obtained from the model Eq. (fT5|l is also shown. The derived model is in 
agreement with the numerical data and the dynamics of the island rotation is recovered. 
We would like to mention here that the time integration for these results has been per- 
formed on a very long time scale compared with the Alfven time. 

Equation (fT5|) shows that the effective island frequency is the result of different 
contributions. In order to investigate the effect of the interchange terms on the island 
rotation. Figure f|T2|) presents the evolution of each frequency for a simulation with 



14 



900 

800 

700 

I 600 
o 

^ 500 
■a 
o 

£ 400 

"O 

c 

f 300 

200 
100- 



- - - numerical position 




^— modei position 




- 








^^^^^ 









0.5 



1.5 

t/T 



2.5 



X 10 



Figure 11. Time evolution of the poloidal position of the island : comparison between 
the model and numerical data for /? = 0.025. 




Figure 12. Time evolution of the frequencies for simulations with /? — 0.025 : [Left] 
Simulation setting k,\~k,i—Q. [Right] Simulation including the interchange terms. 



/3 = 0.025 setting ki = ^2 = (left panel) and for a simulation with /3 = 0.025 
including the curvature/interchange terms (right panel). First, for the two simulations, 
the contributions to the rotation of the equilibrium magnetic field L^^ and of the 
resistivity are weak. Moreover, Figure (fT2|) shows that the frequency dynamics is 
not affected by the curvature terms during the linear regime and during the first quasi 
plateau phase. Actually, during the linear formation of the magnetic island, the rotation 
is controlled mainly by the linear diamagnetic drift while a)* and UexB are weak. During 
the first quasi plateau phase, the nonlinear diamagnetic drift and the E x B poloidal 
flow are strongly generated, and affect the island rotation. The E x B poloidal flow is 
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the most important contribution to the frequency during this regime. However, after 
the first quasi plateau phase, interchange terms affect the dynamics of the frequency. 
When the interchange terms are switched off, the time evolution of the frequencies is in 
agreement with previous results found for a drift tearing mode [1]. After the nonlinear 
generation of the flows, linear and nonlinear diamagnetic drifts cancel each other. As a 
result, the E x B poloidal flow controls the effective frequency of the island, Ut ~ oJexB- 
The right panel of Figure fll2p shows that when interchange terms are included, such 
canceling of the total diamagnetic frequency does not occur anymore after the transition. 
Hence the total diamagnetic drift then provides the main contribution to the island 
rotation. However, after the transition, clearly dut/dt ~ Oojexb/ dt. The asymptotic 
island velocity is enhanced by the curvature terms Ki and K2. 



5.2. Origin of the E x B flow 

Nonlinearly, the E x B flow is generated and affects the rotation of the island. In order 
to investigate its origin, the flow equation (EqU]) is projected on the mode ki for the 
limiting case of ~ d"^. We obtain : 



JLy ^y JLy ^y J Ly 



Hy (16) 



We have defined the E x B poloidal flow as UexB = ^i<9a;0o- So multiplying Eq. 
I6|l by ki and averaging over the current sheet 5, we obtain : 

dt < COExB >5= R{t) + M{t) + V{t) (17) 



where 



R{t) = - ^ [ [ [<P,Vlc^]dydx, 

'JJ-'y J 6 JLy 



^Ly JS JL 



y 



V{t) =^ I I dtHydx 



y 



with R (t) being the Reynolds stress, M (t) being the Maxwell stress and V (t) being 
the viscosity contribution to the E x B flow. 

On the left panel of the Figure (fT3|) . the time evolutions of R, M and V are presented 
for a simulation with /3 = 0.025. As expected, the stresses are nonlinearly generated 
at the beginning of the first quasi plateau phase allowing the growth of the E x B 
poloidal flow. Except at the end of the linear regime where the Reynolds stress is not 
yet generated, the viscosity term is very weak and does not play an important role in 
the generation of the flow. The most important contributions come from the Reynolds 
and Maxwell stresses which balance each other. There is a strong amplification of the 
amplitude of the stresses during the transition. On the right panel of Figure f lT3|) . a 
closeup of the temporal dynamics of R, M and V during the transition time evolutions are 
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Figure 13. Time evolution of M, R and V for /3 = 0.025 ; [Left] Evolution from 
t = Ota to t ~ SOOOOr^i of the stresses over the current sheet. [Right] Zoom during the 
transition of the stresses over the current sheet. 



presented. Note that during this transition, whereas the viscosity term is still weak, the 
evolutions of the Reynolds and Maxwell stresses are complementary. The two stresses 
tend to balance each other during this phase where the flow is crossing the separatrices, 
limiting the level of the generated zonal flow, even if they are both growing in amplitude. 
Once the transition has occured the amplitude of the mean nonlinear brackets in the 
vicinity of the separatrix, R and M, fall. However, in the new dynamical equilibrium 
the resulting E x B poloidal flow persists asymptotically and is driven by the Maxwell 
stress, as a response of the magnetic structure to the new field distribution. 



5. 3. Origin of the nonlinear diamagnetic drift 

In order to investigate the origin of the nonlinear diamagnetic drift, we follow the same 
procedure for the pressure equation (Eq|2]). After projection on the mode ki, we obtain: 

< Cj* >s= dC{t) + dM{t) + D{t) (18) 

where 

dC{t) = -Jf^ I I d4(l),p]dydx, 

y J 6 J Ly 



dM{t) = ^ f I a,[V,Vl^]dyda:, 

y J 6 J Ly 

D{t) = / / dlpdydx. 



SLy 



dC is the contribution to the nonlinear diamagnetic drift of the divergence of the 
convective term, dM is the contribution of the divergence of the Maxwell stress and D 
is the contribution of the diffusivity. 
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Figure 14. Time evolution of dC, dM and D for /3 = 0.025 : [Left] Evolution from 
t = Ota to t — SOOOOr^i over the current sheet. [Right] Zoom during the tansition over 
the current sheet. 



On the left panel of Figure (I14p . the time evolutions of dC, dM and D are presented. 
The three contributions to the nonlinear diamagnetic drift are generated at the beginning 
of the nonlinear regime. During the first quasi plateau phase, dC, dM and D participate 
actively in the generation and growth of the nonlinear diamagnetic drift u*. At the 
transition there is an amplification of the amplitude of the three contributions. The 
right panel of Figure ( I14p presents a closeup of the time evolutions of dC, dM and D 
during the transition. During this transition, the contribution of the diffusivity does not 
grow and is relatively weak compared to dC and dM. Note that both the divergence of 
the convective term and of the Maxwell stress feed the nonlinear diamagnetic drift. We 
observe that in the first phase of the transition, the nonlinear diamagnetic drift is driven 
by dC, the term linked to the advection of the pressure cells. In the following phase 
this term is balanced by the Maxwell divergence stress leading to a stabilization of the 
island dynamics. This shows the importance of the coupling parameter C during the 
transition. In the saturation phase, the origin of comes mainly from the divergence 
of the convective term, dM and D becoming relatively weak. 

5.4- Effect of the /3 parameter on the island poloidal rotation 

In Figure (fT5|) . the effect of the /3 parameter on the rotation frequency is presented. 
The frequencies have been time averaged from t = ISOOOr^ to t = 30000r^. At low /3, 
pressure effects are weak and the situation is close the classical tearing situation. The 
competition between Reynolds and Maxwell stresses produces nonlinearly neither the 
zonal flow nor a diamagnetic drift. The asymptotic rotation frequency Ut increases with 
p. The nonlinear diamagnetic frequency increases almost linearly with /3, but with a 
slope lower than the linear one, allowing a global asymptotic drift of the island in the 
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Figure 15. Effect of (3 on the rotation frequency after the transition. 
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Figure 16. Effete of /? on the stresses after the transition: [Left] Reynolds and Maxwell 
stresses versus /? parameter. [Right] Divergence of the convection and divergence of 
the Maxwell stress versus /3. 



electron diamagnetic direction. Let us note however that the direction of the island 
rotation depends on the value of the viscosity parameter |3]. For any /3, the zonal flow 
contribution to the island drift is weaker than the diamagnetic one. Note also that from 
P ~ 0.015, ojexB decreases, and this is linked to the transition observed in Figure (ITB]) 
where the effect of the /3 parameter on the average stresses is shown (from t = ISOOOr^ 
to t = SOOOOta). The amplification of the stresses leads to the nonlinear generation 
of the mean flows and hence affects considerably the rotation of the magnetic island. 
More precisely, it appears that for low f3 values, the amplitudes of the stresses are very 
weak as in the classical tearing case. This explains why nonlinear diamagnetic drift 
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and nonlinear E x B poloidal flow do not affect strongly the rotation of the island in 
those cases. However, for high /3 regimes, the amplitude of the stresses, in particular 
the Maxwell and convective contribution become more important. At /3 ~ 0.015, a 
transition is observed : flrst, dM ceases to be neglectable compared to dC, second the 
Reynolds contribution R starts to grow, weakening the global E x B flow, as observed 
in Figure (fT5|) . 

6. Summary 

The nonlinear dynamics of a magnetic island in the presence of pressure gradient effects 
has been investigated. This nonlinear dynamics is different from the classical tearing case 
and exhibits a bifurcation. After a linear growth of the island and a flrst quasi plateau 
phase, a transition occurs and the system reaches a new saturated state characterized 
by the flattening of the pressure proflle. We have shown that the dynamics of the island 
during this bifurcation is due to quasi linear effects. The strong generation of a zonal 
flow, due to interchange terms, allows this transition to occur. We have shown that the 
time at which the transition occurs decreases with /3 while it increases with the pressure 
parameter K2. Regarding the poloidal rotation of the magnetic island, a model including 
quasilinear effects has been tested successfully. Before the transition, the rotation of 
the island corresponds to the linear diamagnetic drift. Then, at the transition, the 
rotation is strongly affected by the nonlinear generation of the diamagnetic drift and of 
the E X B flow. We have shown that the asymptotic nonlinear diamagnetic drift is a 
linear function of /3 but does not cancel the linear drift, as previously obtained when 
curvature parameters are neglected. The diamagnetic effect appears to be the dominant 
contribution to the island rotation. We have shown also that the /3 parameter affects the 
magnetic rotation through an amplification of the stresses. We have provided a detailed 
analysis of their impacts on the E x B and diamagnetic drifts for both, the transition 
and the asymptotic regime. At high /3, we find that a Reynolds stress is generated in 
the vicinity of the island and weakens the infiuence of the asymptotic E x B fiow on 
the rotation. 
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